April 20, 2010 0:7 



WSPC - Proceedings Trim Size: 9.75in x 6.5in bilic 



1 



COSMOLOGICAL k-ESSENCE CONDENSATION 



o 

5h _ 

Rudjer Boskovic Institute, 10002 Zagreb, Croatia 
, E-mail:bilic@thphys. irb.hr 



o 



X 



NEVEN BILIC 



GARY B. TUPPER* and RAOUL D. VIOLLIERt 



Centre of Theoretical Physics and Astrophysics, 
University of Cape Town, Rondebosch 7701, South Africa 
* E-maihgary.tupper@uct.ac.za; t E-rnail:raoul.viollier@uct. ac. za 

a -1 

ill 

^q' We consider a model of dark energy/matter unification based on a k-essence type of the- 

ory similar to tachyon condensate models. Using an extension of the general relativistic 
spherical model which incorporates the effects of both pressure and the acoustic hori- 
zon we show that an initially perturbative k-essence fluid evolves into a mixed system 
^ ' containing cold dark matter like gravitational condensate in significant quantities. 

cn 

■^j- ■ The most popular cosmological models such as ACDM model and a quintessence- 

CDM model assume that DM and DE are distinct entities. Another interpretation 
of the observational data is that DM/DE arc different manifestations of a common 
structure. The first definite model of this type was proposed a few years agoP® 
based upon the Chaplygin gas, a perfect fluid obeying the equation of state 



LJ ' A 

P = -~, (1) 
P 

which has been extensively studied for its mathematical properties.^ The general 
class of models, in which a unification of DM and DE is achieved through a single 
entity, is often referred to as quartessenceW^ Among other scenarios of unification 
that have recently been suggested, interesting attempts are based on the so-called k- 
essence^^ a scalar field with noncanonical kinetic terms which was first introduced 
as a model for inflation! 2 -' 

All models that unify DM and DE face the problem of nonvanishing sound speed 
and the well-known Jeans instability. Soon after the appearance of^ andpl it was 
pointed out that the perturbative Chaplygin gas (for early work seep^ and more re- 
cently^) is incompatible with the observed mass power spectrum^ and microwave 
background.^ Essentially, these results are a consequence of a nonvanishing comov- 
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ing acoustic horizon 



dt— . (2) 



The perturbations whose comoving size R is larger than d s grow as S = (p — p) / p ~ 
a. As soon as R < d s , the perturbations undergo damped oscillations. For the 
Chaplygin gas we have d s ~ a 7 / 2 /Hq, where Hq is the present day value of the 
Hubble parameter, reaching Mpc scales already at redshifts of order 10. However, 
as soon as S ~ 1 the linear perturbation theory cannot be trusted. A significant 
fraction of initial density perturbations collapses in gravitationally bound structure 
- the condensate and the system evolves into a two-phase structure - a mixture of 
CDM in the form of condensate and DE in the form of uncondensed gas. 

The simple Chaplygin gas docs not exhaust all the possibilities for quartessence. 
A particular case of k-essence^ is the string-theory inspired tachyon LagrangiarP^I 



C=-V(ip) ^1 - a"" . ( 3 ) 

where 

X = gTw,u ■ (4) 

It may be shown that every tachyon condensate model can be interpreted as a 3+1 
brane moving in a 

4+1 bulkPEH Eq . © is obtained using the stress-energy tensor 
T^ v derived from the Lagrangian with V(ip) replaced by a constant y/A. 

In a recent papei^ we have developed a fully relativistic version of the spherical 
model for studying the evolution of density perturbations even into the fully non- 
linear regime. The formalism is similar in spirit tcE-^and applicable to any k-essence 
model. The key element is an approximate method for treating the effects of pres- 
sure gradients. Here we give a brief description of our method and its application 
to a unifying model based on the Lagrangian Q with a potential of the form 

V(<p) = V n v 2n , (5) 

where n is a positive integer. In the regime where structure formation takes place, 
this model effectively behaves as the variable Chaplygin gas^ with the equation 
of state (JTJ) in which A ~ a 6n . As a result, the much smaller acoustic horizon 
d s ~ a' 7 7 2 + 3 ") i 'Hq enhances condensate formation by two orders of magnitude over 
the simple Chaplygin gas. Hence this type of model may salvage the quartessence 
scenario. 

A minimally coupled k-essence model is described by 



S = I d A x i/—g 



16ttG 



(6) 



where C is the most general Lagrangian, which depends on a single scalar field if 
of dimension m , and on the dimensionless quantity X defined in For X > 
the energy momentum tensor obtained from ([6]) takes the perfect fluid form, 

T^u = 2C X p,p<P,u - Cg^u = {p + p)ujj,u v -pg^, (7) 
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with Cx denoting dC/dX and 4- velocity 



sgn^oA. (8) 



The sign of u M is chosen so uq is positive. The associated hydrodynamic quantities 
are 

p = C(<p,X); p=2XCx(<p,X)-C(<p,X), (9) 

and the speed of sound is defined as^ 

dp 



2 = dp 
8 ~ dp 



s/r. 



dp 



£x (io) 



C x + 2XC X x 



Two general conditions Cx > and Cxx > are required for stability^ and 
causality.^ Now, using ©-(O the <p field equation can be expressed as 



p + m{p + p) + (if - sgn (<p fl )VX)dC/dip = 0. (11) 

Since the 4- velocity (JS) is derived from a potential, the associated rotation ten- 
sor vanishes identically. The Raychaudhuri equation for the velocity congruence 
combined with Einstein's equations and the Euler equation assumes a simple form 

'iU + 3H 2 + a^a^ + 4nG(p + 3p) = ^ ^ , (12) 

V P + P /; M 

where is the shear tensor and = — u^u^ is a projector onto the three- 
space orthogonal to u M . The quantity H is the local Hubble parameter, defined as 
'&% = u v - v . We thus obtain an evolution equation for % sourced by shear, density, 
pressure and pressure gradient. If c s =0, as for dust, Eq. (|12l) and the continuity 
equation comprise the spherical model However, we are not interested in dust, 
since generally c s ^ and the right hand side of (jT2"j) is not necessarily zero. 
In general, the 4-velocity u M can be decomposed as^ 



u^ 1 = (U» + v») I yj\ - V 2 , (13) 

where C/ M = 6ft / •Jgoo is the 4-velocity of fiducial observers at rest, and v^ 1 is space- 
like, with v*- l v^ = —v 2 and U^v^ = 0. In comoving coordinates w M = 0. 

In spherically symmetric spacetime it is convenient to write the metric in the 
form 

ds 2 = N(t, rfdt 2 - b(t, r) 2 (dr 2 + r 2 f(t, r)dn 2 ), (14) 

where N(t,r) is the lapse function, b(t,r) is the local expansion scale, and f(t,r) 
describes the departure from the flat space for which / = 1. We assume that N, a, 
and / are arbitrary functions of t and r which are regular and different from zero 
at r = 0. Then, the local Hubble parameter and the shear are given by 

1 (b l/. \ ^ 2/1 /, \ 2 
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In addition to the spherical symmetry we also require an FRW spatially flat asymp- 
totic geometry, i.e., for r — > oo we demand 

N-H; 6->a(t), (16) 

where a denotes the background expansion scale. 

The righthand side of (fT2"j) is difficult to treat in full generality. As inp^ we apply 
the "local approximation". The density contrast S = (p — p)/p is assumed to be of 
fixed Gaussian shape of comoving size R with time-dependent amplitude, so that 

p(t > r)=p(t)[l + 6 R (t)e- r2 /V R \ (17) 

and the spatial derivatives are evaluated at the origin. This is in keeping with the 
spirit of the spherical model, where each region is treated as independent. Since 
dip = at r = 0, naturally diN — and dib — at r = 0. Hence, 

N(t,r) = N(t,0)(l + O(r 2 )); b(t, r) = b(t, 0)(1 + 0{r 2 )). (18) 

Besides, one finds /q —¥ as r which follows from Einstein's equation G 1 o = 0. 

From now on we denote by H, b, and N the corresponding functions of t and 
r evaluated at r = 0, i.e., H = H(t,0), b = b(t,0) and N = N(t,0). According to 
Q15p. the shear scalar a^a^ vanishes at the origin. Evaluating (|12p at r = yields 
our working approximation to the Raychaudhuri equation. 

We will now apply our formalism to a particular subclass of k-essence unification 
models described by (jSJ). The equation of state is then given by 

P = -« (19) 

P 

and the quantity X may be expressed as 



X(j>,<p) = l-¥t&- = l-<* = l + W . (20) 
P 

The continuity equation, Eq. (|lip. and Eq. (JT5J) evaluated at r = determine the 
evolution of the density contrast. However, this set of equations is not complete as 
it must be supplemented by a similar set of equations for the background quantities 
p and H. The complete set of equations for p, H, ip, 6, p, and H is 

2 

= X(ip,p), (21) 



(If 



^ + 3H(p + p)=0, (22) 

dH 4ttG , , . 

— +H 2 + — (p + 3p) =0, (23) 



(24) 



April 20, 2010 0:7 



WSPC - Proceedings Trim Size: 9.75in x 6.5in bilic 



^ + 3NH(p + p)=Q, (25) 



Ml T. J 



,2 ,4ttG cf (p - p) 



^ + -=- (P + 3p) - 



= 0, (26) 



6 2 i? 2 (p + p) 

where p = p{p, <p) and N = y/X((p,p)/X(tp, p). Eqs. <(STJ> and (22) follow from (JTTJ 
and (|15p. respectively, Eqs. (|lip and (j2"5)l are the continuity equations, and Eqs. 
(|2"3")) and (j2"6")l are the Raychaudhuri equations for the background and the spherical 
inhomogeneity, respectively. 

Now we restrict our attention to the potential (5) . In the high density regime we 
have X ~ 1, and (|2"Tj) can be integrated yielding tp ~ 2/ (3H). Here ~ -ffov^a -3 / 2 
with 17 being the equivalent matter content at high redshift. Hence, V(ip) 2 ~ a 6 ™, 
which leads to a suppression of 10~ 6 of the acoustic horizon at z = 9 for n = 1. 

To proceed we require a value for the constant V n in the potential ([5]). As the 
main purpose of this paper is to investigate the evolution of inhomogeneities we will 
not pursue the exact fitting of the background evolution. Instead, we estimate V n 
as follows. We integrate (|2Tj) approximately with 

X = 1 + w(a) ~ 1 - _ 3 , fi + n A = l, (27) 

as in a ACDM universe^ and we fix the pressure given by ([3]) to equal that of A 
at a = 1. In this way the naive background in our model reproduces the standard 
cosmology from decoupling up to the scales of about a = 0.8 and fits the cosmology 
today only approximately (figure [lja)). 

We solve our differential equations with a starting from the initial adcc = 
l/(z<iec + 1) at decoupling redshift Zd cc = 1089 for a particular comoving size R. 
The initial values for the background are given by 



Pin=Po-3— ; H i n = H o W -g— ; y> in = - — , (28) 



„i ' -iff. 

a dec ,3 - fl » 



and for the initial inhomogeneity we take 



Pin = Pin(l + <5in) , H in = H in (l - , (29) 

where f2 = 0.27 represents the effective dark matter fraction and ^i n = S^a^c) is 
a variable initial density contrast, chosen arbitrarily for a particular R. 

In figure [TJb) the representative case of evolution of two initial perturbations 
starting from decoupling for R = 10 kpc is shown for n — 2. The plots represent two 
distinct regimes: the growing mode or condensation (dashed line) and the damped 
oscillations ( solid line) . In contrast to the linear theory, where for any R the acoustic 
horizon will eventually stop 5r from growing, irrespective of the initial value of 
the perturbation, here we have for an initial <Sfl(adec) above a certain threshold 
S C (R), Sr(o) — > oo at finite a, just as in the dust model. Thus perturbations with 



April 20, 2010 0:7 



WSPC - Proceedings Trim Size: 9.75in x 6.5in bilic 



6 




(a) (b) 

Fig. 1. (a) Evolution of the background in the tachyon spherical model, (b) Evolution of Sr(ci) 
from a doc = 1/1090 for R = 10 kpc, S R (a dac ) =0.004 (solid) and S R (a dcc ) =0.0055 (dashed). 



$R,( a dec) > S C (R) evolve into a nonlinear gravitational condensate that at low z 
behaves as pressureless super-particles. Conversely, for a sufficiently small <5_R(ad ec ), 
the acoustic horizon can stop Sr(o,) from growing. 

The crucial question now is what fraction of the tachyon gas goes into conden- 
sate. IrP3 it was shown that if this fraction was sufficiently large, the CMB and 
the mass power spectrum could be reproduced for the simple Chaplygin gas. To 
answer this question quantitatively, we follow the Press-Schechter procedure^ as 
inP^ Assuming Sn(adec) is given by a Gaussian random field with dispersion o~(R), 
the condensate fraction at a scale R is given by 

where 5 C (R) is the threshold shown in figure^a) . In figureUJa) we also exhibit the 
dispersion 

r°° Ah 

a2{R) ^Ja Y eXP (~ k2R ^ A2{k,adec) ' (31) 
calculated using the variance of the concordance modeP^ 

A* M =co mt (A)V W (_i_)"-\ (32) 

In figure [21(b) we present F(R) for const=7.11x 1CP 9 , the spectral index n s =1.02, 
and the parameterization of Bardeen et alP^l for the transfer function T(k) with 
17^=0.04. The parameters are fixed by fitting (|32|) to the 2dFGRS power spectrum 
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(a) (b) 

Fig. 2. (a) Initial value Snici^cc) versus R for f2 = 0.27 and h = 0.71. The threshold <5 c (ii) is 
shown by the line separating the condensation regime from the damped oscillations regime. The 
solid line gives u(R) calculated using the concordance model, (b) Fraction of the tachyon gas in 
collapsed objects using 8 C (R) and cr(R) depicted in (a). 



dataP^l Our result demonstrates that the collapse fraction is about 70% for n = 2 
for a wide range of the comoving size R and peaks at about 45% for n = 1. 

Albeit encouraging, these preliminary results do not in themselves demonstrate 
that the tachyon with potential (8) constitutes a viable cosmology. Such a step 
requires the inclusion of baryons and comparison with the full cosmological data. 
What has been shown is that it is not correct in an adiabatic model to simply pur- 
sue linear perturbations to the original background: the system evolves nonlinearly 
into a mixed system of gravitational condensate and residual k-essence so that the 
"background" at low z is quite different from the initial one. Because of this one 
needs new computational tools for a meaningful confrontation with the data. 

The tachyon k-essence unification remains to be tested against large-scale struc- 
ture and CMB observations. An encouraging feature of the positive power-law po- 
tential is that it provides for acceleration as a periodic transient phenomenon^ 
which obviates the de Sitter horizon problem.^ 
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